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TECHNICAL NOTE 3318 


ON THE SMALL-DISTURBANCE ITERATION METffi)D FOR THE FLOW 
OF A COSIERESSIBLE FLUID WITH APPLICATION 
TO A PARABOLIC CYLINDER 
^ Carl Kaplan 

SIMMARY 


The assviniptlons of the' Prandtl-Busemann small-dlsturhance method^ 
together wl'th the requirements of continuity and irrotationality, lead 
to a recurs I've system of first-order partial-differential equations. The 
first three sets of equations of this iterati-ve procedure are rewritten 
in con^ilex-'vector form and readily integrated for their particular inte- 
grals . The restilts of the general analysis are then applied to the case 
of subsonic flow past a paraholic cylln^r. This calculation shows that 
the curtailed small-disturbance solution, without the restraining influence 
of a con'fcrol parameter, is unsuitable for the description of subsonic flow 
past the parabolic cylinder. When, howe'ver, the small-disturbance solu- 
tion is developed in powers of the undisturbed stream Mach number Mm as . 

a control parameter and coroared with the solution obtained by means of 
the Janzen -Rayleigh or MM^-esqoansion method, the -two results are identical. 
This agreement shows that the Prandtl-Busemann and Janzen-Raylei^ devel- 
opments are but two different arrangements of the actual solution. Finally, 
the small-distTUbance solution for the parabolic cylinder is examined from 
the point of view of thin-airfoil theory. The series de'velopment of the 
flviid speed at the sm^ace in powers of the ratio of the radius of curvature 
at the •vertex and the abscissa measured from the vertex agrees with the 
resvilts of second-order thin-airfoil theory. Also, a third-order thin- 
airfoil approximation is proposed. 


INTRODUCTION 


The problem of the integration of the equations of conq>ressible flow 
past a prescribed solid boiondary has been treated most often by two approxi- 
mation methods. The first one, initiated by Janzen and Rayleigh, proceeds 
from the incon?)ressible conqilex potential and develops the compressibility 
effects in a series of powers of the undisturbed stream Mach number. It 
is restricted to the subsonic range because the differential equations of 
the process are always of the elliptic type . This method, moreover, is 
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suited particularly for thick bodies and., hence, for relatively small 
critical stream Mach numbers. The second one, the Prandtl-Busemann small- 
disturbance method, proceeds from the undisturbed stream and determines 
the dlstixrbance effects by an expansion in series according to a geometric 
parameter characteristic of the body shape. This method is most stiitable 
for slender bodies for idiich the critical stream Mach numbers are close 
to unity. Very little is known about the limit of convergence of the 
power series employed in these two methods. In a recent paper, however, 
on high subsonic flow past a sinusoidal wall, a plausible argument was 
presented -which indicated that the limit of convergence is coincident -with 
the attainment of local sonic velocity (ref. 1 ) . 

The present paper contains a brief account of Imai's elegant version 
of the Prandtl-Busemann smal 3 .-disturbance method (ref. 2 ) . The original 
Prandtl-Busemann method is based on the following assumption: If e is 

a parameter that characterizes the departure of the profile shape from a 
straight-line segment (for example, thickness, cainber, or an g lR of attack) 
at zero incidence, the -velocity po-tentlEl ^ «-nd the stream function i|r 
can be represented by series in powers of e, the coefficients of -which 
are functions of the flow-plane coordinates x and y. Thus, in non- 
dimensional form: 


9 ! = -X + ejli + e^2 + ^^5^3 + 

^ = -y + 6^1 + e^2 




(1) 


where the undisturbed stream is directed from right to 3 eft. On the basis 
of this assumption, there foUo-ws from the general second-order nonlinear 
compressible-flow equation for ^ or by means of a conparlson of coef- 
ficients, a recucrsive system of second-order differential equations for 
the coefficients of the individual powers of e . The first one is a Laplace 
type of equation and the ones that follow are of the Poisson type, the 
right-hand sides of -which are composed of pre-viously determined functions. 


Imai’s version of the Prandtl-Busemann method proceeds from the set 
of first-order differential equations for ^ and if that resuilts from 
the requirements of continuity and irrotationality. Thus, 





(2) 
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Then, hy means of equations (l) and the following expression for — , 



a comparison of the various powers of e yields a recursive system of 
first-order equations for • • • similar to 

equations ( 2 ) . ■ These pairs of equations can be expressed in complex- 
vector form and readily integrated for their particular integrals. The 
symbols tised in the preceding equations are defined as follows; 

x,y rectangular Cartesian coordinates in flow plane 

^ velocity potential 

i|r stream function 

u,v velocity conponents in direction of x- and y-axes, respectively 
p density of fluid 

p„ density of fluid in undlstxirbed stream 

q fluid speed 

Mgg Mach mmOaer of undisturbed flow 

7 ratio of specific heats at constant pressure and volume 

The quantities x, y, if, u, v, and q are all nondimensional and the 
subscripts denote differentiation with respect to the desig na ted variable. 

Whereas in the Janzen-Baylelgh method the expansion of the complex 
potential in powers of is always possible, in the PTandtl-Busemann 

method the e^anslon of i|r (or 0 ) in the form of equations (l) cannot 
be guaranteed a priori. Indeed, experience has shown that only in cases 
where the prescribed profile has no stagnation points so that the assump- 
tions of the small-dlstTirbance method are adhered to strictly can the 
es^ansions for ^ and ijr, as indicated in equations (l), be valid (refs. 5 
and . Most profiles of aeronautical interest, however, have rounded 
leading edges and therefore possess stagnation points, in the neighborhood 
of -vrtilch the asstmiptlon of small deviation from undisturbed flow is clearly 
violated. In such cases, sooner or later, terms of the form e^log € 
must appear (ref . 5) • 5 ?he procedure then is as follows : The velocity 
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potential and stream functions are still assumed in the form of equa- 
tions (1), hut with the parameter e now considered as a dumriQr syiabol 
which serves only to regulate the course of the Iteration process . VJhen 
the stream function only is considered^ the approximation func- 
tions iJjjCn ^ 2) are assumed to satisfy the boundary condition that they 

vanish on the contour, are regiilsr in the region of flow outside the bound- 
ary, and their derivatives vanish to a sufficient degree at infinity. The 
first pair of terms in the expansion usually represents the Prandtl-Glauert 
linearized approximation. It is further assumed that, except for a Bm'n 1 1 
region in the nei^borhood of stagnation points, the function and its 

derivatives are small compared with and its derivatives for ail m < n. 
The question whether this assimiption is satisfied can be answered only after 
the have been calculated. 

As borne out affirmatively in the only two cases thus far calculated, 
namely, the elliptic, cylinder and the circular-arc profile (refs . 5 an^ 8 ) , 
it seems reasonable to conjecture that, if it were possible to calculate i|fn 
to any order and .if each term were developed in powers of M»^, the formal 
arrangement in series of these powers would yield precisely the Janzen- 
Raylel^ solution for the same profile. The example chosen to illustrate 
the general analysis of the present paper and to verify the preceding 
conjecture is the flow at zero incidence past a parabolic cylinder . This 
profile is especially amenable to treatment by means of both the small - 
disturbance method and the Janzen-Raylei^ method. Moreover, because the 
parabolic cylinder represents a magnified picture of a stagnation region, 
it is particularly well suited for a critical examination of the utility 
of the smal 1 -disturbance method for the calculation of the flow in the 
neighborhood of stagnation points . It does not follow, even though the 
preceding conjectvire is verified, that the small-disturbance method provides 
a siiltable approximation in the neighborhood of a stagnation point. 


AMALYSIS 


The development of the nondlmensional fluid speed q. is obtained 
with the aid of equations (l) and ( 2 ) . Thus, 

q 2 - 1 = -24 ~ "^2^ + 

2^^(9^1x^2x ^^ly?^2y " ?^5x) + * * * 
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The corresponding development for ^ is obtained by means of equa- 
tions (3) and that is. 


^ = 1 - + I 


(1 + - -4. 


2x 









+ 


(5) 


A comparison of coefficients of individual povers of e in equations (2) 
leads to the following sets of recursion formulas for the first three 
steps : 


P^ljc = 

“ “^ix 


( 6 ) 


^^2x = ^2y - I ' ^Ix^ly j' 


M7) 


J^2y = -^2x + 
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3 ^ 3 x = ■ 4 ' 3 y - 


^2y ■ I “ 


I - 1(1 + 


lx ’^'ly 


I ^^ly^’l'ly ■*■ ^2x’l'iy + ^2s^2y 


hv = ’V ■*■ 


- if + ’^^js'nx^'l'a. 


1 2 

- 01y iflx + ^2x^'lx 


By the introduction of the con^lex notation 

X + iPy = z 


X - iPy = z 


9^n + J tn = ’Wh 

9 ^n - 7 fe = % 
P 


> (8) 
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and with the aid of the symbolic relations 


9 oc Sz ^z" 

Sy ySz bz j 

equations (6), (7), and (8) become, respectively. 



( 9 ) 


'' 2 z “ “ 5 


(<,-1) 


p _ 2 

+ 2ctw. W-_ + crw,_ 
Iz Iz Iz Iz 


( 10 ) 


-jr = i *-2.)- - - 2 z *- 2 .*- 2 s)* 


i| 

k 


2 - p2(4cr - 1 ) 


- ''a/) + ^ - ■"as) + 


5 8 J (’'a^ ''izf 


(iJ.) 


p -imj-2 _ _ y+1 Mdo 

where = 1 - Mw and <r = 1 + 

4 p 2 


Note that the right-hand 


sides of equations (lO) and (U) are composed of, respectively, doiible 
and triple products of previously determined perturbation quantities. 
Equations (9), (lO)^ and (U) are first-order complex-vector equations 
with Vjj as dependent variable and z and z as independent variables . 

They can be integrated in a straightforward manner. Thus, equation ( 9 ) 
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means that is a function of z only. Then, the general integral of 

equation (iD) is the following egression: 


VA — — — H33 


(CT - 1 )ZW- 


^ 2 
Iz 


+ 2 owiWi^ 




+ Hz) 


( 12 ) 


■vdiere F(z) is an arbitrary function of z to be determined by the 
boundary conditions. 

From equation (l 2 ), the general Integral of equation (U) is as 
follows; 


''3 = 4“^^ 




“'ll + 


i(a<r - 1) [a (3 + p2) - + 

(or - 1 ) [1(3 + p 2 ) - - 1) [20(3 + - 

(5 + 342)]&j^ 5 + p2) . 1] + 

I o(ff - + I 

* 

2^1^12^12 + WiTTj/ ) ' - I “ l)^lzPz + 

o(^wiF 2 + Wi2.F^+ cLzj + G(z) 


(13) 
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-vHnere G(z) is an arbitrary fimction of z to be determined by the 
boundary conditions . 


It is noteworthy to remark that in the original version of the small- 
disturbance method the differential equations 'corresponding to equa- 
tions (9)^ ( 10 ), and (U) are of the second-order Laplace and Poisson 
types. Therefore, the arbitrary functions added to the particular 
integrals can be functions of either z or "z only. Past experience, 
however, has shown that certain terms in the third approximation give 
rise to singularities in the region of flow (ref . 6) . These spurious 
singularities must be con^jensated for by the addition of suitable func- 
tions of z and of 'z in compliance with the boundary conditions at the 
sxu’face and at infinity. In the present version of > the small-disturbance 
method, the arbitrary functions added to the particular integrals can 
only be functions of z. QJherefore, because both functions of z and 
functions of z are necessary for the removal of external singularities 
in the presence of a solid boundary, the complication of apparent singular- 
ities in the external region cannot occur in the present circumstances. 

As a final remark, note that the form of the particular integrals obtained 

for subsonic flow remains the same for supersonic flow with = Moo^ - 1 


and cr = 1 


7 - 

p2 ‘ 


SUBSONIC PLOW PAST A PABABOLIC CYLINnER 


Before proceeding to the calcixlation of the subsonic flow, the 
Incon^ressible flow past a parabolic cylinder will be derived. Thus, 
consider the transformation 




•vrtiere z = x + iy and 5 = 5 + Then the- real and imaginary parts’ 
give 


X = 



( 14 ) 


y = 241 


> 
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and hence, hy the elimination of the variable ti. 


(15) 

Thus, the CTurves for | constant are parabolas vhose foci are at the 
origin (fig. 1 (a)) . In order to obtain the flov past a parabolic 
cylinder | = |q, vith undisturbed velocity from right to left, it is 

observed that the nondimensional stream function is given by 


if = -2(1 - io)n 

The con^)lex potential w = ^ + iilr is then given by 

V = -(^2 _ 

and the complex velocity by 


dv 

dz 


u - iv 



Thus, 




u = -1 + 


I + Tl 


V = 

■ + tj2 


(16) 


n 
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Now, from eq\iations (l 4 ) , 


52 = 


X + 


\/x2 + y2 



-X + y/ x2 + y2 
2 


> 


The polar equation of the parabola is 


(17) 


r 



1 + cos 0 


vhere 


X = r cos 0 


and 


y = r sin 0 


Hence, from the second of equations 


(17), on the parabola 


I = lo. 


Ti = ± lo tan I 

Therefore, from equations (16), 


u 


-sin^ — 


and 


T = 


± — sin 0 
2 
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or 


q. = \[\^ + = sin ^ 

According to Bernoulli's theorem, the pressure coefficient Cp is 


Cp - 



cos 


2 6 
2 


Figure 2 shovs the -velocity and pressure-coefficient dls-bributlons 
along the upper svirface of the paraholic cylinder. Note the monotonlc 
character of q. and that the undisturbed s-tream speed is reached at 
x=y = - oo(e = jt). This heha-vior makes -fche paraholic cylinder a partic- 
ularly good shape for the examination of the Prandtl-Busemann small- 
distTxrhance me-thod in -the nei^horhood of a stagnation point. 

Consider now the soibsonlc case; let 


X + ipy = 



(18) 


or 



y = 2 |ti J 


(19) 


Elimination of 


Tj' from these equations yields 
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p P 

Thus, the curves for | constant are parabolas with foci at x = > 

o *** 

y = 0 with the focal distance independent of the stream Mach num- 
ber Mm (fig. 1 (b)). Note, further, that equation ( 20 ) results from 
equation (18) with the coefficient of iT] completely arbitrary. The 


choice of was made in order that tj be identical with the ti of the 

p 

incon^jressible case at the surface. Thtis, from equations 03 ), 


|2 = 


X + \/x^ + 


2 p^ 




( 21 ) 



-X + ^x2 + p2j^ 
2 


j 


Also, on the parabola I = §0 (^iS* ^(l^)),- 


X = -M„,^|q^ + r cos 0 


and 


y = r sin 0 


where 



1 + cos 0 


Introducing these expressions into the second of equations ( 21 ) yields 


n 


= ± Iq tan I 


as in the inconipresslble case. 
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Now, the first approximation ^ is a function 

of z = X + ipy and therefore of S = I + ^ n only. Conrparison with 

3 

the incoaipressihle case shows inmifidiately that 


^1 = 9^1 + I = 2lo5 


Thus, 


= -2(1 - |o)n 


(22a) 


which satisfies the boundary condition that ilf = 0 for I = lo that 

at infinity the disturbance velocities vanish for points not near the 
parabola. The esqiresslon for the velocity potential becomes 






+ 26 oI 


(22b) 


Now, let the positive sense of describing the parabolic boundary be 
counterclockwise with the positive neormal direction inward. Then, if ds 
and dn represent elements of arc and inward normal, respectively, the 
egressions for the normal velocity in terms of the stream function and 
the velocity potential are as follows (fig. 1 (b)) ; 


^00 

qn = --^ 


( 23 a) 


2\/lo + 


On - “ 


2(32|o2 


— ^ ■ ■ (lo^ + - W»^lon5^Ti (23b) 

+ — — 
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Clearly, from equations (22a) and (23a), the normal velocity vanishes 



flnri ( 23 b) . Thus, in general, the hotindary conditions cannot he satis- 
fied simultaneoTisly for both ^ and \|r. In the present version of the 
small-disturbance method, the function ^ does not represent the velocity 
potential of the flow but is tcfcilized only for purposes of notation and 
ease of calculation of the stream function ilr. 

How, with wi = 2go5 and z = eq.uation (l2) yields 


+ F z(0 


The arbitrary function F is detemdned from the boundary condition 
that •>lf 2 = 0 on the parabola 5 = “5 + 2 |q. Thus, 


1 t 2 


F ' - 

(ff - 4* Ikr Y log ^ 



where the right-hand side of this equation is a pure ima gi n ary. Also, 

I.P. -1.P. £ 

f £2 

I.P. i= -i.p. i 
I.p. log f = -I.p. log ^ 


I.P.F =» -I.P.F 
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Hence, on the hotindary. 



1 ) 


~ _ ka ^ ~ - 2ff log S 


The expression on the right-hand side involves the varlahle ^ only and 
is regular thm ughout the field of flow^ therefore. 


TP_ 1 . 2 M„ 

^-4 So ^ 


(a . i) (S - 2 {q) _ _ ^ S - 2to . 20 log E 


and 


1 . 2 M« 

’^ 2 =“r^o ~2 


(CT - 1) 


12 - (S - 2|o)^ ^ 


S' 


S 


Thus, 


^2 = I.P.pW2 = Mc^lo^CS - lo)^ 


g 4 - 1 
+ Ti^ 


+ 2(g - 1) 


P^S(l + So) 


2\2 


+ 11 ^) 


J(24) 


This expression for '^2 satisfies the boundary condition that it vanishes 
at the surface of the parabolic cylinder I = Iq and alo n g the x-axis 
(y « Tfj ss 0) outside the boundary. In a dd it ion, it vanishes at infinity. 

The expression for the third approximation is obtained from 

equation (15) - In the following,, there are listed some of the individual 
terms of equation (ij) and the correspo n d ing terms of the arbitrary func- 
tion g|z( 53 chosen in such a way that it is regular everywhere in the 
field of flow and that 1^5 vanishes along the parabola I = So 
ainng the x-axis y w q = 0 ; 
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^Iz^^lzz 


lo^ tS - atg)' 


2^'^ 55 


2p» j: 




zw^— ^ 

p4 p 




21^1 


■dz = ± log 5 


Sk!^aog-L 

^ - 26o 2|o 


v,v, 2^ 2|o3^ 

1 pl^ ^2 

zw 5 = lo^ L 
P^ 53 


2ap3 s - 2lo 

P^ 52 

U - 2|,)2 

P^ 53 


2So^ 1 


/'wi^3a^^_£io£ 1 

^ p^ 5 p4 5 


M 


Wi2‘ 


2\ ^ . 2|o 3 |3 ^ ^ (s - 2 So)3 


P^ 






Izz 


2gp^l^ 

p^ 53 


P^ 

Io£Lz_2|o 

p^ 

2|o^ (g - 2 |o)2 


P'^ 


> ( 25 ) 


2V. w w -lio£l 

‘^l^lz^lz ^ = “7- 7 

p^ 5 p4 5 


- P 2^0^ S 

= — r 

p4 p 


2Sq^ S - 2 |q 
p^ 52 



CBIH 
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Thiis, 




-a(ka - 1) - 


- 1) [ia(5 + p2) . (5 ^ 


- 2U 


(2CT - i)[a(3 + p2) . log X + I 

-•S - 2|o ® 2io ^ 




- 1 ) + 


2(a - 1) j^(3 + p2) _ 


5 - 2|, 


2 - 


|{o - 1) [a(3 + 92) - l]| 1 . yi^%(a . .) ■» - ^toA _ 


So3*^<!(a-l) 


(g - 2gp)g 


5^ 

S - 2|o 


- 2ff 


- DSo® 11- - + 2(a + 1)1^ (S - 26o)^ 


Ho - 1)|„2 L^_?k + 2 („ 


5' 


+ 0 


-(0 - i) X ^^o , + inj 1 + 

£2 £ 


-— ~ — log X 
5 - 26o 2So 


+ 20 


2 " 
- 1)^ + (o + I)i2 + 2 
? 52 J 


I 
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This expression for G is regular throughout the field of flow, since 

the apparent pole at ^ = 2 |q is canceled hy the zero of log — ^ there. 

26o 

Then the expression for w^ is given by 


|lfe,2ff(a - i)ll±^Llil2L + |l(ff _ 1)[20(5 + p2) _ (5 + 33a)] _ 


(Aio^i 


, f - (£ - 2 |o )2 


+ 4(ff - i)[ff(3 + + 


I ff(3ff - l)K»^ " 3 + P^) - ^ - 3o%4.^| ^ + 

^ ■ i (t ““ i - rAi: ifc)) * I ^ ■ 

a(tr + i)|a + . (5 _ 2io)^ + arj^(ff - + a(a + i)i| + 

^[l+(S-2io)| -ff(o 


Le + (s - 2io)| - «r|(ff - 1) ^^^ ~ - £ (S “ 2gp) _ ^ a - 26 q) - £ 


^ log -i- + — log -i 

5 2|o s - 26o 26 o 


^ - a>|^o - 


•tSiI 

d L’ vt 

’ £7 r n 1' 

vj 



20 


NACA TN 5518 


and 


- I + 5lo|2 - 2I0H - Wlo5) - 

pVE®(ii£5 + x55^e 2 - aoio^s - aoio3) + ^'*(i5|5 + ooioi^ - ioio®E - - 

il^(5E + So)l + I - " 

-* + r) '- 


ap^l(>H^ + toE - *»So®) + n*^(6E + Eo)j + £ 0 ^® — -— - jo^‘ 


(»2|2 + r,^y 


|(ff - 1) ^7 + p 2 ) - (5 + 3 P^ij| + 3 lo) - (31 + Eo)n^ - 

+ 1) As z l ° ^ L-(3gg|a - na) - 1 - D |[»(n + 5pS) -^t- 

(p2e2 + , 2)3 I 


0M,2eoI + |Eo^u 4(® - 1)(T + P®) + ^ /g - g - 

J(p2|2 + ,,2)2 2 L J(Bat2 + 


E - Eo 


(p2t2 + r^) p(| _ a|„)a + ^ 


n log 


p2|2 + i)a 

4 P^Eo" 


P2(|2-2S„QH.,2 ^^_^ H 

PEo PE 


(a6) 


This expression for satisfies the houndary condition that it vanishes 

at the sxirface of the parabolic cylinder I = |q and along the x-axis 

(y w T) = 0 ) outside the houndary. The congolete e:gpresslon for the stream 
function ijr, inclusive of the third approximation, is then given hy the 
sum of equations (22a), ( 24 ), and (26). 

For the purpose of conparison vith the Jans»n-Raylei^ solution (to 
he derived in the following section), the expression for the fluid speed 
at the houndary will now he obtained in the form of a power-series devel- 
opment in Moo^- Symbolically, 




NACA TN 5518 


21 


and 


8y 


J: r, A + p2i A 

2 . St) 


^{f-v + n‘ 


Also, in terms of the relevant quantities. 


^=l-^ikly + ^ + t|,^2) - i|r2y + t|T2x>lr2x + 


P2 


^^0y^2y - ^^3^ 


+ • • • 


Then, if only the terms that involve the Mach number up to ere 

retained. 


= -2(g - lo)tl + 2Mooh/(t - lo)' — - ■- + 

|2 + if- 


loT log 


|2'-t- t ]2 


(|2 + q2 _ 2|ol)tan“^ 1 


&o^(l “ lo) 


(1^ + T]^) (1^ + Tl^ - Ii-IqI + 


and 


£a=l-_ 




P + Tl' 

Uov, along the surface of the cylinder I = |q^ “lost of the interest 

centers. 


2 ^ „2 


1 =- 


2p(pho^ + 11^) 


<^|)p= 


l=lo 
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JAHZEET-RAYTiBIGEg MmOD FOR PARABOLIC CYLmDER 


The solution of the problem of siibsonic flow past a paralolic cylinder 
ly means of the Janzen-Raylei^ method follows along the same lines as that 
ly means of the smaU-disturhance method. Thus^ 
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and 




t = llfo + + - • • 


> 


(28) 


where and ilfQ are, respectively, the incompresslhle velocity 

potential and stream function. Again, from equation (3)4 expanding in 
powers of yields 



(29) 


Then, from equations ( 2 ), a conqjarison of the various powers of 
yields the following recursive system of first-order equations for 
9^0^ ’I'o^ ^1^ l^l^ • • •* 


Now, let 



^^ly = 


2 ^ox ^ " 


(0, 


'ox 



X + iy = z 
X - iy = ¥ 

^n + ifc = ■«h 

^n “ 


(30) 


(31) 
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Equations (30) and (31) then hecome 


= 0 


(52) 


and 


1 — 1 — 2 
"Iz ]j^ OZ JL oz oz 


( 35 ) 


Equation (32) niftans that Vq Is a function of z only, and the general 
Integral of equation ( 33 ) is 


= “7 Wq + f Wq2 / Vq^ 6S + H(z) ( 3 ^) 

4 4 J 

Tdiere H(z) Is an arhitrary function to be determined according to the 
boundary conditions. 

For the parabolic cylinder, 

z = 5 ^ (£ = I + n) 


and 


Wo = -(C - lo)^ 


From equation (^) then 


-2 - 5 

^1 = I ^oI - j lo^ I + J So Y - lo^ I + I ^ log I + h[z(£)] 



4F 
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The arbitrary funcblon H is determined, from the boundary condition 
that = 0 on the parabola ^ = -^ + 2 |q. Thus, 


2i>lri = i lo(5 - 5) - ^ lo^Clog 5 - log 0 + i I 






I lo^f| 2og S - i log n + H - H 


where the ri^t-hand. side of this equation Is a pure imaginary. Also, 
I.P.f(5) = -I.P.f(5); hence, on the boundary, 


I.P.H = I.P. 


i to5 - 1 lo® ^ - i So - 


2S, 


So® ,®^° - i log J. 

C 2 S - 2So 21, 


The expression on the ri^t-hand side involves the variable ^ only and. is 
regular throughout the field of flow; therefore. 


■2 =■ I So5 


i I 2 iQg _£ 1 I " 2lo)2 

„ So , Sq 


2 lc 


t 


io2 - i -iol. log ^ 


2 5 - 21 , 


21 . 


5 
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and 


»1 = I 5 o « + t ) 


-2 

1 . 5 




- (5 - 2 | o )' 




2 ^ ^ - 3h . + 1 1^54 log J i — log X) (35) 

° C 2 ° k 2|o ^ - 2|o 2i ' 


Then, 


- 60S -.1 to^ ^ ^ - 

^60 


|2 + tj 2 


lo ( 1 ^ - 2|ol - Ti^)log + 2 ti ||2 + tj 2 - 2|ol 

^lo 


and 


2 | o 2 ) tan-l a 


(1^ + Ti^) (|2 + tj 2 > l^i^i + 41^2^ 


( 56 ) 


T | r ^ = - - . Jp . l ^. + 


+ Ti ^ 


^ o 2 

lo^ log 5 


2 | ol ) tan”l a 


lo^Cl - lo ) 


( 1 ^ + + f ?- 


(37) 
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Usually, most of the Interest lies In the velocity distribution on the 
houndaiy. Thus, substitute I = Iq Into equation (36) and. in the 
-E.P .(5 - then at the boundary 






2 21 

5 + 


log 




+ 2 |o tan' 


-1 JL 




The "velocity of the fluid, is gi-ven by 


a = 



Hence, 


q = 




<2n - uj- 


2 


* 2 2 


Moo' 


2 ^o' 




fo1 + (lo^ - -3- 

1,1 2 So 


This equation is in agreement vith equation (27) . Thus, the small- 
disturbance solution vhen expanded in powers of yields precisely 

the Janzen-Eaylelgh result. Now, at the upper surface, T = Iq 
therefore. 
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and 



cos^ ^ 
2 


1 + f 

4 


1 + 3 sin^ + sin 0 


6 cos 6 - 



(59) 


Table I lists -values of q. and Cp for Moo = O.5 o-ver the upper sur- 
face of the cylinder and figure 2 shows the corresponding graphs. 


Since the completion of this paper the attention of the author bas 
been drawn -bo a recent calculation by Tmal of the Janzen-Baylei^ solution 
for a parabola including terms in Moo^ (ref. 7) • ^ error in sign in 

reference 7 has been corrected. (Note last term in equation (38)) . 


DISCUSSION OF ANALISIS 


The main concern of this paper has been the presentation of Imai's 
elegant -version of the smal 1 -disturbance method and its application to 
the problem of -two-dimensional compressible flow past a parabolic cylinder. 
The example of the parabolic cylinder -weus chosen for the diial pui5>ose of 
illustrating the results of -fche general analysis and for comparison -with 

the Moa^ -expans ion or Janzen-Raylei^ method of solution. 

One of the basic assumptions of the small-disturbance method is that, 
except for a small region in the neighborhood of the nose, the deri-va- 
ti-ves of T|r^ must be less than those of where m < n. A numerical 

comparison at the surface of -the parabolic cylinder shows, howe-ver, 
that ' > >t ^py o-ver a large portion of the surface . Presijmably, 'if it 

■were possible to calculate to any order, the condition that ilfnx < ifaix 

for n > m -would be satisfied from some definite -value of m on-ward . 
Althou^ -the solution as obtained in the present paper satisfies the 
bOTjndary conditions at the surface and at Infinity for each step of the 
iteration process, -violation of the foregoing basic assun5>tion renders 
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the curtailed solution useless for the calculation of the velocity and 
pressure distributions at the surface. Nevertheless, some" interesting 
results vith regard to thin-alrfoil theory can be obtained from the 
expressions for 1^2^ and obtained from equations (22a), , 

anfl (26) . Thus, by foUoving the ideas of Van Dyke (ref. 8), the velocity 
at the surface of the parabolic cylinder is developed in a series of 


powers of 



the ratio of the radius of curvature at the vertex and 


the abscissa measured from the vertex. For this purpose, consider the 
following expression for the fluid speed q in terms of the deriva- 
tives 


q = 1 + ^ " ’*^2y " %y 2 ‘*‘2 ’’'ly^ 2 ’*'2y 

♦ix’t'ac + Vay) * ^ ®*3y’''2y 




where only those terms have been retained which involve the ratio — £ 

inclusive of the second power. From equations (22a), (2k), and (26), 
the required expressions for the derivatives of ijfg, and i|Tj at 

the surface of the parabolic cylinder are as follows: 
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u _ ^°h r 2 lo- 


i|rpx = “ Moo^C^^ + 

^2 ^3 




+2y=|“«^ (‘^ + l)^-*-2(cr-3)p2^ 


^3y = M»2 i o2(i . 5p2) _ |(5 „ p2) + ^1 . 5p2) 


iji(ff - 1)(T + P^) + ^ 


' I 1 -11^ 

2 log 2^ ^ 

M I E 

. P^ P^ 

\ n T , 


Siibstitution of these escpressions Into equation (!«)) then yields 
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q. = 



1 + -L(iwM„2 _ I'N 

1 + (ff - l)MU^fl + p2 + ) 

if \ ^ an" 


\ / 

\ / 




i +1)2-1 5^|2 (j 2(1 - 5p2) + 4ffp2 + g - ZJ^ + 

p2 3 p2 L J 


2 p2 



- l)(7 + p2) 





(42) 


The first two terms on the right-hand side of this equation agree with 
the results of second-order thln-alrfollr theory (ref. 8 ). The third 
term presumably wo old he obtained from a third-order thln-alrfoil theory. 
Although this term has been written in a form to suggest that it is of 


order 



It is of interest to note that the esspression 



— 1 

tan 



can be expanded as follows; 
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The -third-order -thin-alrfoil result can then he -written as f oUowb : 


q = 1 - ± 


hA 

1 + (ff - l)M„^ 

W 

- 


s/at 

52 


Mco^jff(d - 1)(7 + 



m 


The small circles in figure 2 designate points calculated by means of 
this equation -with Mo, = 0-5* The first circle corresponds -to 0 = 120° 

2 


for -which -the parame-ber 
^2 


2 E, 


equals — . If -the terms that in-volve 

3 




are included^ the magnitude of q. becomes greater than unity 


e-ven for -values of 0 corresponding to points relatively far from the 
nose. For example, q. = 1.0132’ at the point corresponding to 0 = 120°. 
These spurious -values of q. indicate that -fihe fourth approximation of 
the small-disturbance method must con-tribute additional terms of the 
.2 


order 


(¥) 


Langley Aeronautical Laboratory, 

National Ad-vlsory Committee for Aeronautics, 
Langley Field, Va., Sep-tember 21, 195^* 
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TABLE I 

■ 7 EL 0 CITY AMD PEIESSIIEIE-TOEEFICIENT DISTEIBOTIONS 


■ 

00 

= 0 

= 

D.50 

a 


9 


0 

0 

1.0000 

0 

1.0625 

5 

.01^36 

.9981 

.0403 

1.0607 

15 

.1305 

.9830 

.1207 

1.0459 

30 

.2588 

.9330 

.2402 

.9971 

45 

.3827 

.8536 

.3573 

.9185 

60 

.5000 

.7500 

.4705 

. 814-7 

75 

.6088 

.629k 

.5781 

.6915 

90 

.7071 

.5000 

.6782 

.5565 

105 


.3706 

.7687 

.4184 

120 

.8660 

.2500 

. 84-72 

.2866 

135 

.9239 

.1465 

.9115 

.1707 

150 

.9659 

.0670 

.9596 

•0795 

165 

. 992 h 

.0170 

.9896 

.0207 

180 

1.0000 

0 

1.0000 

0 
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